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Abstract 

We consider the branching ratio of b — ► 57 in gauge mediated supersymmetry 

breaking theories. Useful bounds on the parameter space of these models are 

derived from the experimental bounds on b — > sj. Constraints on masses of 

NLSP are presented as a function of tan/3 and M/A for fi < and fi > 0. 
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There has been a tremendous recent interest in the phenomenological implications of 
gauge mediated supersymmetry breaking theories [1|-|5| • These theories are characterized by 
gravitino as the lightest supersymmetric particle(LSP) and have signatures for supersymme- 
try which are distinct from the usual minimal supersymmetric standard model(MSSM). The 
interest in these theories is especially heightened because of the possible explanation of the 
peculiar event seen at CDF with final state containing e + e~77 |§ and missing transverse 
energy. A successful explanation of this event requires the neutralino to be the next to 
lightest supersymmetric particle (NLSP). Whether or not this explanation withstands the 
test of time, it would seem important to examine in detail the mass constraints on NLSP 
that ensue in these models from the rare decay b — > 57 where SUSY contributions occur in 
one loop diagrams. Although some preliminary work exists in special cases H, ours is the 
first comprehensive study in the full allowed parameter space. 

We shall analyze b — > 57 in the gauge mediated SUSY breaking model in the whole range 
of combinations of A and M (where A is the effective SUSY breaking scale and M is the 
messenger scale), the number n of (n 5 + n 5 ) pairs, and a range of tan/3 starting from small 
values to large values. We take A ~ O(100TeV) since soft SUSY breaking scalar masses 
are then of the order of the weak scale. The parameter is M > A. There could be a large 
hierarchy M » A ||, however the upper bound of the gravitino mass ~ IKeV restricts 
the M/A < 10 4 @. The number of pairs n is restricted to 4 by the perturbative unification 
of the gauge coupling. The 10 representation can be included by noting that one n w + n w 
pair corresponds to n=3. 

The induced gaugino and scalar masses at the scale M are ||: 

*dO-.(^)2gQA. (i, 

m\M)=2(n)f(±) gfc jCj (!*M) 2 A2 (2 ) 

where ki and Cj are 1,1,3/5 and 4/3, 3/4, and Y 2 for SU(3), SU(2), and U(l), respectively. 
The values of Cj apply only to the fundamental representations of SU(3) and SU(2) and are 
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zero for the gauge singlets. 

We use the exact messenger-scale threshold functions [|[ 
1 + x 

d( x ) = — t-\og(l + x) + (x->-x), (3) 
x z 



, s. 1 + ar 
f( x ) 



x' 2 



log(l + x) - 2Li 2 ( — — ) + -Li 2 ( 



— . + (,T — > -x) , (4) 

l + xj 2 \l + xj\ v ; v ; 

rather than the limiting values, /(0) = g(0) = 1 and /(l) = 0.7 and g(l) = 1.4. We shall 
require that electroweak symmetry be radiatively broken. We use ot s = 1.20, sm 2 9 w = 0.2321 
and a = 1/127.9 at the weak scale as the gauge coupling inputs. We first go up to the 
messenger scale M with gauge and Yukawa couplings, and fix the sparticle masses with the 
boundary conditions ([!]) and (|2]). We next go down with the 6x 6 mass matrices for the 
squarks and sleptons to find the sparticle spectrum for large as well as small tan (5. We use 
the RGE's given in the reference || . Also, we do not choose a particular model for fi. Note 
that the soft Higgs mass parameters m 2 H and m 2 H2 from eq. (0), along with the ratio of the 
vacuum expectation values tan/3(= Vi/v 2 ) uniquely specifies where /i is the coefficient 
of the bilinear Higgs term in the superpotential.We will use two extreme values of tan/3 
equal to 3 and 42 (for n = 1) for illustration. It is interesting that in the case of tan (3 = 3 
neutralino is the NLSP(next to lightest SUSY particle) for all values of M > A; however 
for tan/? = 42, either the stau or the neutralino is the NLSP depending on whether \i and 
M/A are small or large. For n — 2, however stau can be NLSP even for the small tan (3 for 
the lower ratios of M/A. When stau is the NLSP, these models can not explain "the" event 
in Tevatron. However the scenarios can have other interesting collider signatures both in 
Tevatron as well as LEP. Calculation of h — >• 37 amplitude involves the coefficients of short 
distance photonic and gluonic operators c^{M w ) and c$(M w ). Effects of QCD corrections 
to two loops is then carried out. For the Standard Model these calculations are given in 
Ref. ||10|| . Recent attempts to make partial corrections in the next to leading order |TTJ do 



not indicate any large terms. Contributions from various supersymmetric contributions are 
given in a generic form in Ref. [ . We use our calculated mass spectrum and the couplings 



to calculate b — >• 37 rate. The results depend on A, M, tan/?, n and sign of //. We shall 



use more physical variables tan (3, /x, sign of /x, M/A and n. Our figures are for the minimal 
case of n = 1, and we shall remark on the situation for higher n. The total amplitude has 
contribution from the W-loop, charged Higgs (H*) loop, chargino {x^) loop, neutralino 
loop and the gluino (g) loop. We find that the neutralino and the gluino contributions to 
the amplitude are less than 1% in the whole range of parameter space. The charged Higgs 
contribution adds constructively to the W-loop contribution. The chargino contribution 
can occur with either sign, but is generally much smaller than the Higgs contribution. An 
exception is when tan (3 is large and fi < and chargino interference opens up the allowed 
parameter space. 

We consider branching ratio vs /x for either sign of /x in two different scenarios: (a) 
tan/? = 42 and M = 1.1A and M = 10 4 A; (b) tan/? = 3 and M = 1.1A and M = 10 4 A for 
n=l. The two different values of M form the boundaries of the envelope of parameter space 
that would be traced by any relation between A and M. We exclude the case M = A since it 
produces a massless scalar in the messenger sector. We have used the particle data values for 
the CKM matrix elements and also have imposed the constraint | V^V^ 2 / (V^ 2 = 0.95±0.04 

In Figure la corresponding to scenario (a) we display b — > 57 branching ratio as a 
function of |/x| for /x > 0. Solid lines represent M = 1.1 A and the dashed lines M = 10 4 A. 
CLEO bound 1 x 10" 4 < Br(b -> 37) < 4.2 x KT 4 at 95% CL clearly rules out the smaller 
values of /x. We find \x > 720. When the branching ratio is 4.2 x 10 -4 the NLSP (stau in the 
case of M = 1.1 A) mass is 182 GeV. Fig lb displays the scenario(a) for /x < 0. The chargino 
destructive interference for larger values of M/A does not yield any useful constraint on the 
parameter space. In this parameter space, lighter mass NLSP solutions correspond to stau 
as NLSP. The extreme left ends of the curves correspond to a bound on the lightest slepton 
mass ~ 65 GeV. 

In Fig 2a, 2b we consider two cases in scenario (b), for positive and negative /x respec- 
tively. The left ends of the curves correspond to the lowest chargino mass bound which we 
have taken to be around 75 GeV. In these cases, chargino contribution is very small. The 
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variation with M is also small. We then find that the CLEO constraint leads to /x > 460 
GeV for positive /x and |/x| > 426 GeV for negative /x when M = 1.1 A. 

Constraints on /x can be translated into bounds for the masses of supersymmetric par- 
ticles. We are interested in displaying these bounds for the masses of NLSP. In Fig 3a. we 
display the lower bound on neutralino mass as a function of tan (3 for ji > and /x < for 
two limiting values of M, i.e. M = 1.1A and M = 10 4 A. Again as tan/5 becomes larger, 
chargino interference in /x < case for M = 10 4 A removes any useful bounds. Though the 
variation of the ratio of M/A does not produce much difference in the branching ratio for 
the lower values of tan /3, but it has bigger effects in the bounds for the NLSP. In Fig 3b we 
display the lower bounds on NLSP when tan/? > 31. Four cases considered are similar to 
Fig 3a. However, we do not have any CLEO bound on the NLSP when /x < and M = 10 4 A 
throughout the range of tan (5 displayed. The solid curve and the dashed curve corresponds 
to the bounds on stau mass which is the NLSP and the dot-dashed curve corresponds to the 
bounds on the neutralino mass, which is the NLSP in this case. The bound on /x, for the 
positive values of /x is a monotonic function increasing from 460 GeV for tan/3 = 3 to 720 
GeV for tan/3 = 42 when M = 1.1 A. For M = 10 4 A, the bound increases from 480 GeV to 
920 GeV in the same range of tan (3. These large values of /x raise a problem of fine tuning. 

For n — 2, the constraint on the NLSP mass is higher, e.g. for tan/3 = 3, lowest mass for 
the stau (since it is the NLSP) allowed by CLEO data would be 128 GeV when /x < and 
M = 1.1A, however for M = 10 4 A the lowest mass allowed for the NLSP is 72 GeV and it is 
neutralino. Almost the same thing happens for /x > 0, the lowest stau mass allowed is 136 
GeV for M = 1.1 A and the lowest neutralino mass is 115 GeV for M = 10 4 A. We also have 
the same conclusion for the tan (3 = 42, i.e the lowest mass for the stau (since it is the NLSP) 
allowed by CLEO data would be 83 GeV when /x < and M = 1.1 A, however for M = 10 4 A, 
there is no bound on the NLSP (stau), for /x > 0, lowest stau mass(NLSP) allowed is 183 
GeV for M = 1.1A and the lowest stau mass (NLSP) is 201 GeV for M = 10 4 A. We have 
assumed the value of 175 GeV for the top running mass. The results are rather insensitive 
to this mass. A variation of 5 % in mass results in the change of branching ratio of less than 
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1%. 

In conclusion, we have used the CLEO bound on the branching ratio for b — > 57 to limit 
the parameter space of the gauge mediated supersymmetry breaking models. We have found 
useful bounds on masses of NLSP. We also have found that for positive //, irrespective of 
tan (3, jj, is restricted to large values. Since this raises the problem of fine tuning, our analysis 
shows that gauge mediated model generally favors negative \i solutions. When \i is negative 
the available parameter space increases with the ratio of M/A. In the near future, with an 
improved bound on the branching ratio for b — > 57, it will be possible to put more severe 
constraints on the parameter space. 



This work was supported by Department of Energy Grant DE-FG03-96ER-40969. 
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Fig. 1 : Plots for b — > 57 branching ratio as a function of for tan (3 = 42. 
Solid lines correspond to M — 1.1 A. 
dashed lines correspond to M — 10 4 A. 

a) for pi > 

b) for fx < 0. 

Fig. 2 : Plots for 6 — > S7 branching ratio as a function of for tan (3 = 3. 
Solid lines correspond to M — 1.1A. 
dashed lines correspond to M — 10 4 A. 

a) for /j, > 

b) for /i < 0. 

Fig. 3 : 3a) Bound on NLSP (neutralino) mass as a function of tan (3. 

Curves from the top (around tan (3=5) correspond to cases: (i) /i > and M = 1.1A, 
(ii) fj. < and M = 1.1A, (iii) // > and M = 10 4 A, (iv) /i < and M = 10 4 A. 
3b) Bound on NLSP (stau for the solid and dashed lines and neutralino for the dot- 
dashed line) mass as a function of higher range of tan (3. 

Curves from the top (around tan/5=32) correspond to cases (i) /i > and M = 1.1A, 
(ii) fi > and M = 10 4 A, (iii) /i < and M = 1.1 A. 
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